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We demonstrate that near the threshold, the tttt scattering amplitude contains a pole with the 
quantum numbers of the vacuum - commonly referred to as the a - and determine its mass and 
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relies on a straightforward calculation based on the Roy equation for the isoscalar S-wave. 
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According to the Particle Data Group, the lowest res¬ 
onance in the spectrum of QGD carries angular momen¬ 
tum £ = 0 and isospin 1 = 0. The state is listed as 
/o(600) and is usually called the a. It manifests itself 
as a pole on the second sheet of the isoscalar 5'-wave of 
TTTT scattering. We denote this partial wave amplitude by 
tg(s). The numbers for the pole position found in the lit¬ 
erature cover a very broad range. For recent reviews, we 
refer to Him In fact, since such a state is not easily 
accommodated in the multiplets expected for qq bound 
states and glueballs, some authors question its existence. 

All of the pole determinations we are aware of rely on 
models and, moreover, use specific parametrizations to 
perform the analytic continuation. In the present paper, 
we instead rely on an equation which has been shown to 
follow from first principles, the dispersive representation 
of the partial wave amplitude to(s) due to Roy 0: 


to complex values of s and use this extension to (a) prove 
the existence of a second sheet pole close to the threshold 
and (b) determine the position of this pole within rather 
small uncertainties. For this purpose, we only need the 
particular equation quoted above. The explicit expres¬ 
sion for the kernels occurring therein reads: 
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It amounts to a twice subtracted dispersion relation. 
Crossing symmetry implies that both subtraction con¬ 
stants can be expressed in terms of the S'-wave scattering 
lengths: a = a^, h = {2 — 5ao)/(12M^). The integral 

describes the curvature generated by the S- and P-waves 
below A and the so-called driving term do(s) collects the 
dispersion integrals over the higher partial waves {£> 2), 
as well as the high energy end of the integral over the S- 
and P-waves. 

Similar equations hold for all other partial waves. 
Those for the S- and P-waves amount to a set of cou¬ 
pled integral equations, which strongly constrain the low 
energy properties of these waves 0 Q ■ Previous work 
on the Roy equations concerned the behaviour on the 
real axis. In Ref. 0, a crude estimate of the mass and 
width of the a was given, but as emphasized in several 
reviews (see e.g. 012), this estimate relies on a para¬ 
metric extrapolation off the real axis and is thus subject 
to a sizeable systematic uncertainty. 

The present paper closes this gap. We show that the 
domain of validity of the Roy equations can be extended 


The first term in Kq{s,s') accounts for the contribu¬ 
tions generated by the right hand cut. The remainder of 
Kq{s, s') describes those contributions from the left hand 
cut that are also due to Imtg, while those from lmt\ (P- 
wave) and Imtg (exotic S'-wave) are booked separately. 

For our analysis, it is essential that the dispersion in¬ 
tegral is dominated by the contributions from the low 
energy region: because the Roy equations involve two 
subtractions, the kernels Kn{s,s') fall off in proportion 
to 1/s'^ for large s'. Note that the contributions from 
the left hand cut play an important role here: Dropping 
these, Ko{s,s') reduces to the first term, which falls off 
only with the first power of s'. Taken by itself, the con¬ 
tribution from the right hand cut is therefore sensitive 
to the poorly known high energy behaviour of ImtQ(s'), 
but taken together with the one from the left hand cut, 
the high energy tails cancel. 

Since the decomposition into partial waves is useful 
only at low energies, the dispersion integral over the S- 
and P-waves in Eq. m is cut off, the high energy tail 
being included in ^[((s). As discussed in Ref. a, d 0 (s) 
is dominated by the contribution from the P-wave res¬ 
onance 72 ( 1270 ). Using the narrow width approxima¬ 
tion, expanding the relevant kernel in inverse powers of 
s' = and retaining only the leading term of order 
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FIG. 1: Behaviour of the amplitude near threshold. 
1/s'^, we obtain 

^ 5(»-4Ma(lU + 4Af;)r,.^„ 

The experimental values for mass and width are Mf^ = 
1275.4 ± 1.2 MeV, = 158.5 ± 4.4 MeV 0. 

We have performed a detailed evaluation of the driving 
term, which exploits the Roy equations for the D-, F- 
and G-waves and accounts for the contributions from the 
high energy tail. The calculation shows that, in the low 
energy region we are considering in the present paper, 
the contributions from partial waves with £ > 2 as well 
as those from energies above 1.4 GeV are negligibly small 
and the narrow width approximation works: for A > 1.4 
GeV, the above simple formula is within the uncertainties 
attached to the full result, which affect the outcome for 
the pole position only by 1 or 2 MeV. 

The most important feature in the low energy region is 
the occurrence of an Adler zero. To leading order of chi¬ 
ral perturbation theoryMbe amplitude is given by Wein¬ 
berg’s formula of 1966 0, 

t°(s) = (2s-MM327rF2) . (4) 

In this approximation, the zero occurs at s = \M'^. FiglU 
shows the behaviour of the amplitude - calculated from 
Eq. ~ near the threshold and explains why the theo¬ 
retical predictions are so precise in this region: by far the 
most important contribution stems from the subtraction 
term, which is linear in s and is fixed by the two scat¬ 
tering lengths Oq, a§. Accordingly, the values of Og and 
Gq represent the most important ingredient of our calcu¬ 
lation. Since theory predicts these very accurately 0, 

aj] = 0.220 ±0.005 , ag =-0.0444 ± 0.0010 , (5) 

the uncertainties in the subtraction term are very small. 
The main experiments concerning the scattering lengths 



FIG. 2: Domain of validity of the Roy equations. 


IISII1II3 are all in good agreement with these predic¬ 
tions, but provide a stringent test only for the first one. 
On the other hand, as will be discussed below, the experi¬ 
mental information does suffice to evaluate the dispersion 
integrals in Eq. m to good accuracy. 

The Roy equations represent the partial wave projec¬ 
tions of the fixed-t dispersion relations obeyed by the 
TTTT scattering amplitude. These relations are valid if the 
variable t is contained in a Lehmann-Martin ellipse with 
foci at t = 0 and t = — s' and right extremity at 

t = r(s'), where s' > AM^ is the variable of integration. 
In the Mandelstam representation, the size of the ellipse 
is limited by the singularities due to the double spec¬ 
tral function. For AM^ < s' < 20 M^, the correspond¬ 
ing expression for r(s') reads r(s') = 16s'/(s' — AM^), 
while for s' > 20 M^, the right extremity is at r(s') = 
As'/is' -16 M^). 

The partial wave projection involves the values of 
the amplitude on the interval t = ^ — s)(l — z), 

6 < z < 1. For the Roy equations to be valid at the 
point s, the corresponding interval must be contained in 
the intersection of the ellipses characterized above. The 
boundary of the domain G where this is the case is shown 
as a solid line in fig. [3 using pion mass units. The bound 
for r js' ) established on the basis of axiomatic field theory 
[lilLll4| is only slightly weaker. The dash-dotted line in¬ 
dicates the corresponding domain of validity of the Roy 
equations. On the real axis, this domain reduces to the 
range —AM^ < s < 60 obtained in Ref. 0. 

The values of the 5'-matrix element 

5°(s) = 1 - 2x/4M2/s-l (s) (6) 

on the second sheet can be calculated from those on the 
first sheet: unitarity implies the relation |l 5|| 

s^isY^ = 1/s^oYY ■ ( 7 ) 

The Roy equation thus automatically also specifies the 
function Sq(s) on the second sheet, in the same domain 
of the s-plane . In particular, the amplitude contains 
a pole on the second sheet if and only if Sq(s) has a zero 
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on the physical sheet. So, all we need to do is numerically 
evaluate Eq. o for complex values of s in the domain G 
where it has been shown to hold and find out whether or 
not Sq{s) has zeros there - standard routines that solve 
equations numerically immediately provide the answer. 

As discussed in Ref. 0, the Roy equations for the 
S and P-waves provide very firm control over the low 
energy behaviour of the imaginary parts occurring in 
Eq.(^. The numbers given in the following are based 
on a new analysis of these equations which will be de¬ 
scribed elsewhere. We could just as well have used the 
results in Ref. 0 - the numbers would barely change. 
For our central representation of the scattering ampli¬ 
tude, the function S'g(s) has two pairs of zeros in G, one 
at s = (6.2±i 12.3) M^, the other at s = (51.4±z 1.4) M^. 
These are indicated in fig. [3 which may also be viewed 
as a picture of the second sheet - the dots then repre¬ 
sent poles rather than zeros. In the following, we work 
with the complex mass Too- = ^ Tc, defined as the 

value of ■’/s at the pole on the lower half of the sec¬ 
ond sheet. For the central solution, the pole occurs at 
TOct = 441 — i 272 MeV, not far from the place where the 
tj was resurrected ten years ago Ea 

The second pole represents the well-established res¬ 
onance /o(980). In order to study the behaviour of 
the amplitude there, we have extended the analysis of 
Refs. 0 0 to higher energies. The BES data on the de¬ 
cay J/V' —*■ (/)7r7r El have clarified the structure in this 
region. In this reference, the sharp drop in the elasticity 
r;0o(s) at the KK threshold is parametrized by means of 
the so-called Flatte formula, which describes the interfer¬ 
ence between the KK cut and the pole from the /o(980). 
Given the elasticity, the Roy equations determine the 
phase shift. It turns out that the solution closely fol¬ 
lows the input: the phase shift differs from the phase of 
the Flatte formula only by a slowly varying background. 
It does not come as a surprise, therefore, that our calcu¬ 
lation confirms the position of the pole in and we do 
not discuss this further. 

Finally, we estimate the uncertainty to be attached to 
the result obtained from our central representation of the 
scattering amplitude. The Roy equations imply that if 
the two subtraction constants Og and Og are treated as 
known, the low energy properties of the isoscalar S'-wave 
depend almost exclusively on a single parameter, which 
may be identified with the value of the phase at 800 MeV 
and which we denote by 6 a = 5g (800 MeV). Accordingly, 
it is convenient to represent the position of the pole as 

mrj = mo + mi Aqq + m 2 Aal + m 3 ASa , (8) 

Aa)i = (aj] - 0.22)/0.005 , Aag = (ag -k 0.0444)/O.OOl, 
A(5^ = ((5 a -82.3)/3.4. 

The term toq represents the value of the complex mass 
obtained if the scattering lengths are fixed at the central 
values of the prediction 0 , while the phase at 800 MeV 
is set equal to the central value obtained from the phe¬ 
nomenology of the phase difference between and 6 ° 0. 


The coefhcients mi, m 2 and m 3 describe the sensitivity 
of the result to a change in these variables - for the small 
changes of interest, the linear approximation is adequate. 

The residual noise in toq is small, for the following 
reasons, (a) The integral over Imtg is dominated by the 
region below the KK threshold. If Og, a§ and 6 a are held 
fixed, the Roy equations constrain the integrand very 
strongly there, (b) The contribution from the P-wave 
is dominated by the p. Since the experimental informa¬ 
tion about this state is excellent, the integral over Im t\ is 
known very well, (c) In Im tg, the uncertainties are larger, 
but the entire contribution from this wave is small, (d) 
As stated above, the uncertainties in the higher partial 
waves and in the contributions from the high energy tail 
of the dispersion integrals affect the result for mu only 
at the level of 1 or 2 MeV. In our opinion, the estimate 
Too = (441 ± 4) — i (272 ± 6) MeV generously covers the 
uncertainties. For the other coefficients in Eq. 0, we ob¬ 
tain mi = —2.4 -k i 3.8, m 2 = 0.8 — i 4.0, m 3 = 5.3 -k i 3.3 
(values in MeV). 

The result for the pole position does not change signif¬ 
icantly if the imaginary parts are evaluated with a phe¬ 
nomenological representation of the data. We illustrate 
this with the parametrization of the scattering amplitude 
proposed in appendix A of Ref. [T^. Evaluating the dis¬ 
persion integrals over Im tg, Im t\ and Im tg with the cen¬ 
tral parametrization of that reference, we find that the 
pole occurs at 445 — i 241 MeV. Although this represen¬ 
tation differs significantly from our central solution, the 
result for the pole position agrees with formula 0 : for 
the values Og = 0.23, Og = —0.048 and 5a = 90.9° that 
correspond to this parametrization, the formula yields 
Too- = 447 — i 242 MeV. 

This confirms that the position of the pole from the 
tj is indeed controlled by three observables. Only one 
of these, Og, is known experimentally to good precision. 
For the second one, Og, the sharp theoretical prediction 
in Eq. 0 was recently confirmed by an evaluation on the 
lattice: Ug = —0.0426(6)(3)(18), where the three brackets 
give the statistical, systematic and theoretical errors, re¬ 
spectively (cf. also 113). To stay on the conservative 

side with the third parameter, we use 5a = 82.3° -40 . 
Compared to the uncertainty from this source, the noise 
in the term mg is negligible. With the theoretical predic¬ 
tion for the scattering lengths, formula 0 then leads to 
our final result 

= 441 tg® MeV, T,, = 544 Izf MeV . (9) 

We conclude that the same theoretical framework that 
leads to incredibly sharp predictions for the threshold pa¬ 
rameters of TTTT scattering 0 also requires the occurrence 
of a pole on the lower half of the second sheet, with the 
quantum numbers of the vacuum, not much above the 
threshold, but quite far from the real axis: the width of 
the a is larger than the width of the p by a factor of 
3.7. The parametric extrapolation used in Ref. 0 led to 
somewhat higher values, for the mass as well as for the 
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width - the difference amounts to about 1 standard de¬ 
viation. The uncertainty in the result stems almost 
exclusively from 5a- The range adopted for this observ¬ 
able covers all of the phenomenological parametrizations 
we are aware of. (Energy-independent analyses have a 
broader scatter of values, but even the most extreme 
is less than 2 a away from our central value.) It can 
be reduced substantially if the data that underly these 
parametrizations are compared with the solutions of the 
Roy equations. We intend to describe this elsewhere. 

The experimental information concerning Uq is meagre, 
but since the real part of the coefficient m 2 is very small, 
this does not affect the value of Mg-: as far as the real 
part of the pole position is concerned, the result remains 
practically unchanged if the theoretical predictions for 
the scattering lengths are replaced by the experimental 
constraints on these [ij, US • 

Many of the determinations of the mass and width of 
the a listed by the Particle Data Group neglect the left 
hand cut. In the language of Eq. m, this approximation 
amounts to replacing the kernel Kq by the first term in 
Eq. lO and dropping Ki, K 2 as well as dp. Eor defi¬ 
niteness, we fix the subtraction constants a and b such 
that the “exact” and approximate representations agree 
at the threshold and at --/sa = 800 MeV. At the ener¬ 
gies of interest, the difference is then well described by 
the parabola ~ c(s — 4:M^){sa — s), with c = 

0.5GeV“"'. Removing this term from has a rather 
drastic effect: the pole then occurs at — 500 — i 260 
MeV. The amplitude obtained in this manner cannot be 
taken at face value, because it violates unitarity: drop¬ 
ping the left hand cut necessarily also distorts the imag¬ 


inary part. We do not pursue this further. The above 
expression for the corresponding curvature shows that 
the left hand cut cannot be neglected - the pole is not 
sufficiently far away from it for this approximation to be 
meaningful. 

It is difficult to understand the properties of the low¬ 
est resonances in terms of the degrees of freedom of the 
quarks and gluons. The physics of the a is governed by 
the dynamics of the Goldstone bosons: The pro perties of 
the interaction among two pions are relevant . A qual¬ 
itative emlanation for the occurrence of the a was given 
in Ref. |3, on the basis of current algebra, spontaneous 
symmetry breakdown and unitarity. The properties of 
the resonance /o(980) are also governed by Goldstone 
boson dynamics - two kaons in that case. It would be of 
considerable interest to apply the above analysis to the 
Roy-Steiner equation for the Ktt S'-wave with I = \- 
This should clarify the situation with the k. 


Acknowledgments 

We are indebted to David Bugg for numerous com¬ 
ments concerning the issues discussed above, in particu¬ 
lar for providing us with detailed information about the 
structure of the amplitude near the KK threshold. Also, 
we wish to thank Peter Minkowski and Wolfgang Ochs 
for informative discussions. This work is supported by 
Schweizerischer Nationalfonds, by the Romanian MEdG 
under Contract CEEX 05-D11-49 and by the EU “Eu- 
ridice” program under code HPRN-CT2002-00311. 


[1] W. Ochs, in Hadron Spectroscopy: Tenth International 
Conferenee, edited by E. Klempt, H. Koch and H. Orth, 
AIP Conf. Proc. No. 717 (AIP, New York, 2004), p. 295. 

[2] M. R. Pennington, arXiv:hep-ph/0509265 

[3] D. V. Bugg, in Hadron Spectroscopy: Eleventh Inter¬ 
national Conference on Hadron Speetroscopy, edited by 
A. Reis, C. Gobel, J. de Sa Borges and J. Magnin, AIP 
Conf. Proc. No. 814 (AIP, New York, 2006), p. 78. 

[4] S. M. Roy, Phys. Lett. B 36 (1971) 353. 

[5] J. L. Basdevant, C. D. Froggatt and J. L. Petersen, Nucl. 
Phys. B 72 (1974) 413. 

[6] B. Ananthanarayan, G. Colangelo, J. Gasser and 
H. Leutwyler, Phys. Rept. 353 (2001) 207. 

[7] G. Colangelo, J. Gasser and H. Leutwyler, Nucl. Phys. B 
603 (2001) 125. 

[8] S. Eidelman et al. [Particle Data Group], Phys. Lett. B 
592 (2004) 1. 

[9] S. Weinberg, Phys. Rev. Lett. 17 (1966) 616. 

[10] S. Pislak et al. [BNL-E865 Collaboration], Phys. Rev. 
Lett. 87 (2001) 221801 Phys. Rev. D 67 (2003) 072004. 

[11] B. Adeva et al. [DIRAC Collaboration], Phys. Lett. B 
619 (2005) 50. 

[12] J. R. Batley et al. [NA48/2 Collaboration], Phys. Lett. B 
633 (2006) 173. 

[13] A. Martin, Scattering Theory: Unitarity, Analyticity and 


Crossing, Lecture Notes in Physics, Vol. 3, (Springer- 
Verlag, Berlin, 1969). 

[14] G. Mahoux, S. M. Roy and G. Wanders, Nucl. Phys. B 
70 (1974) 297. 

[15] For a derivation, see for instance Ref. 

[16] We thank Bachir Moussallam for this remark. 

[17] N. A. Tornqvist and M. Roos, Phys. Rev. Lett. 76 (1996) 
1575. 

[18] M. Ablikim et al. [BES Collaboration], Phys. Lett. B 607 
(2005) 243. 

[19] J. R. Pelaez and F. J. Yndurain, Phys. Rev. D 71 (2005) 
074016. 

[20] S. R. Beane, P. F. Bedaque, K. Orginos and M. J. Savage 
[NPLQCD Collaboration], hep-lat/0506013 

[21] T. Yamazaki et al. [CP-PACS Collaboration], Phys. Rev. 
D 70 (2004) 074513. 

[22] R. Kaminski, L. Lesniak and K. Rybicki, Eur. Phys. J. 
directC 4 (2002) 1. 

[23] V. E. Markushin and M. P. Locher, in Workshop in 
Hadron Speetrocopy, Frascati Physics Series, Vol. XV 
(Laboratori Nazionali di Frascati, Rome, 1999), p. 229. 

[24] P. Buettiker, S. Descotes-Genon and B. Moussallam, Eur. 
Phys. J. C 33 (2004) 409. 



